function res = funcPhi(kh,dd)

res = besselh(0, kh*dd);



